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ABSTRACT 

We present a detail analysis on the spectral lags of the short gamma-ray bursts (GRBs) 
and compare them with that of the long GRBs by using the CGRO/BATSE GRB 
Catalog. Our sample includes 308 short GRBs and 1008 long GRBs. The light curves 
of these GRBs are in 64 ms time bin and they have at least one long and intense pulse, 
which satisfies ST > 0.512 seconds at c = la and c max > 6cr, where ST is the pulse 
duration, c is the photon counts, and er is the standard error of the background. We 
calculate the cross correlation function (CCF) for the light curves in 25-55 keV and 
110-300 keV bands and derive the spectral lag by fitting the CCF with a Gaussian 
model. Our results are as follows. (l)The spectral lag distribution of the short GRBs is 
significantly different from that of the long GRBs. Excluding the statistical fluctuation 
effect, a proportion of ~ 17% of the short GRBs have a negative spectral lag, i.e., 
the hard photons being lag behind the soft photons. We do not find any peculiar 
features from their light curves to distinguish these bursts from those with a positive 
spectral lag. We argue that a more physical mechanism dominated the hard lag may 
be hid behind the morphological features of the light curves. This should be a great 
challenge to the current GRB models. We notice that this proportion is consistent 
with the proportion of short GRBs correlated with nearby galaxies newly discovered 
by Tanvir et ah, although it is unclear if these short GRBs are indeed associated with 
the sources originated at low redshift. (2)While the spectral lags of the long GRBs 
are strongly correlated with the pulse durations, they are not for the short GRBs. 
However, the ratios of the spectral lag to the pulse duration for the short and long 
GRBs are normally distribution at 0.023 and 0.046, respectively, with the sample 
width, indicating that the curvature effect alone could not explain the difference of 
the spectral lags between the two types of GRBs. The hydrodynamic timescales of the 
outflows and the radiative processes at work in GRBs might also play an important 
role as suggested by Daigne and Mochkovitch. 
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1 INTRODUCTION 

The phenomenon of gamma-ray bursts (GRBs) is still a 
great puzzle in astrophysics, although much progress has 
been made (see reviews by Zhang & Mezaros 2004; Piran 
2005). Both the light curves and the spectra among GRBs 
are enormously various. There is competing evidence that 
multiple GRB populations exist. The large sample detected 
by CGRO/BATSE was identified as two types of GRBs, 
long-soft and short-hard GRBs, separated by a duration of 2 
seconds (Kouveliotou et al. 1993). It is well established that 
the long GRBs are associated with energetic core-collapse 
supernovae(e.g., Galama et al. 1998; MacFadyen & Woosley 
1999; Bloom et al. 1999; Hjorth et al. 2003; Thomsen et 


al. 2004). Although it is well known that the short GRBs 
are phenomenologically different from the long ones, such 
as variability time scale (Liang et al. 2002a; Nakar & Pi¬ 
ran 2002), fluence (Liang et al. 2002b), pulse interval and 
number of pulses in a burst (Nakar & Piran 2002), hardness 
ratio (Kouveliotou et al. 1993; Qin et al. 2000, 2001; Dong 
& Qin 2005; Qin & Dong 2005), and spectrum (Ghirlanda 
et al. 2004; Cui, Liang, & Lu 2005), the issue that whether 
or not their energetics and host galaxies are different from 
the long GRBs remains unknown being due to poor local¬ 
ization and lack of afterglow observations before Swift era. 
Some groups paid much effort to search for the afterglows 
for the short GRBs, but no detections were made (Kehoe 
et al. 2001; Gorosabel et al. 2002; Hurley et al. 2002; Klotz, 
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Boer & Atteia 2003). Swift satellite (Gehrels et al. 2004) 
now is ushering us in a new era in short GRB study, and 
rapid progress was made in the first Swift operation year 
(Gehrels et al. 2005; Hjorth et al. 2005a, b; Berger et al. 
2005; Fox et al. 2005). The precise localization of short GRB 
050509B by Swift (Gehrels et al. 2005) led to the first red- 
shift measurement and host galaxy connection (Bloom et al. 
2005) for this class of GRBs. Among four short GRBs with 
known-redshift so far, GRBs 050509B (z=0.2248, Bloom et 
al. 2005), 050724 (z=0.257; Berger et al. 2005), and 050813 
(z=0.72; Prochaska et al. 2005) coincide with the early-type 
stellar population with no or little current star formation, 
favoring mergers of compact object binaries as the progen¬ 
itors of the short GRBs. one case, GRB 050709 (z=0.16), 
occurs at a late-type dwarf galaxy with a star formation 
rate exceeding 0.5 Mq. yr -1 (Prochaska et al. 2005). These 
breakthroughes rapidly improve our understanding on the 
nature the short GRBs. 

In this paper we present a detail analysis on the spectral 
lags of short GRBs. The observed spectral lag is a common 
feature in high energy astrophysics (e.g., Norris et al. 2000; 
Zhang et al. 2002). It is found that soft lag, soft photons 
being lag behind the hard photons, is dominated in long 
GRBs (Norris et al. 2000; Wu & Fenimore 2000; Chen et al. 
2005; Norris et al. 2005). The lag is also strongly correlated 
with the luminosity (Norris et al. 2000) and the jet break 
times in afterglow light curves (Salmonson et al. 2002). We 
here present a study on the spectral lags of the short GRBs 
and compare them with that of the long bursts by using 
CGRO/BATSE GRB Catalog. The data and our GRB sam¬ 
ple selection are described in §2. We present our analysis 
method in §3, and compare the spectral lags in Short and 
Long GRBs in §4. We find that the spectral lags of the long 
GRBs are correlated with the pulse duration (ST), while the 
spectral lags of the short GRBs are not. We study the re¬ 
lation between the spectral lag and the pulse duration and 
the ratio between these two quantities in §5. It is interest¬ 
ing that about one-third of the short GRBs have a negative 
lag (hard lag, hard photons being lag behind the soft pho¬ 
tons). We present a further analysis on these short GRBs in 
§6. Discussion and conclusions are presented in §7 and §8, 
respectively. 


2 DATA AND SAMPLE SELECTION 

The sensitive GRB survey made by CRGO/BATSE presents 
a large and homogenous 1 samples for both long and short 
GRBs. The light curves of these bursts are in four energy 
bands, i. e., 25-55 keV, 55-110 keV, 110-300 keV, and > 300 
keV. We use the light curve data in a time bin of 64 ms, 
which are concatenations of three standard BATSE data 
types, DISCLA, PREB, and DISCSC 2 . We take the DIS- 
CLA data before GRB trigger as the background of each 
GRB. We fit the DISCLA data by a linear function, and 
then subtract the background with this linear model. Our 
long and short GRB sample include 1008 and 308 GRBs, re¬ 
spectively. The light curves of these GRBs have at least one 

1 It is recently suggested that the BATSE short GRB sample 
may have two components(Tanvir et al. 2005). 

2 ftp://cossc.gsfc.nasa.gov/compton/data/batse/ 


long and intense pulse, which satisfies ST > 0.512 seconds 
at c = la and c ma x > 6<r, where c is the photon counts and 
a is the standard error of the background. 


3 ANALYSIS METHOD 


The cross correlation function (CCF) has been widely used 
to measure the time lag of two light curves in two energy 
bands (Link, Epstein, & Priedhorsky 1993; Cheng et al. 
1995; Band et al. 1997; Norris et al. 2000; Li et al. 2004; 
Chen et al. 2005; Norris et al. 2005). Assuming that {vi} 
and {v 2 } are the light curves in two different energy bands, 
with N data points in a time bin of At, the CCF as a func¬ 
tion of r = kAt is defined as 


CCF(t = kAt-,vi,v 2 ) 


EfT i(f)n 2 (t + r) 
NsJ<j' vl a' V2 


(1) 


where 



i= 1 


( 2 ) 


The CCF as a function of r for two light curves with the 
same profile but having a lag of to symmetrically peaks at to. 
The observed light curves are mixed with noises. The CCF 
hence may not symmetrically peak at to. We do not simply 
read off the lag from the peak of the CCF curve. In order to 
reduce the scattering caused by noise , we use a Gaussian 
function to fit the CCF curve, and take the peak of the 
Gaussian as the lag. The lag derived by this way should be 
more precise than that directly read off from the CCF. It can 
be a fraction of the time bin. A positive to indicates a soft 
lag, corresponding to an earlier arrival time of higher energy 
gamma-ray photons than that of lower energy photons. 


4 COMPARISONS OF SPECTRAL LAGS IN 
SHORT AND LONG GRBS 

We calculate the lag between the light curves in 25-55 keV 
and 110-300 keV bands. The lag derived by the CCF is an 
average time lag (t) for the whole time domain. It is found 
that some GRBs, such as GRB 960530 and GRB 980125, 
have significantly different spectral lags in early and late 
epoches(Hakkila & Giblin 2004). In principle, one could not 
measure multiple lag components for two light curves by 
the CCF. Therefore, we also calculate the spectral lag of 
each pulse in the light curves. The comparisons of the dis¬ 
tributions of the spectral lag measured from the whole light 
curves (t), spectral lag of the first pulse (ti th), and the max¬ 
imum spectral lag 3 among pulses within a GRB ( T max ) for 
the long and short GRBs are shown in Figure 1. We summa¬ 
rize the features of these distributions in Table 1. The null 
hypothesis that these distributions for the long and short 
GRBs are drawn from the same parent is tested by the K- 

5 test, and the results are also tabulated in Table 1. From 
Figure 1 and Table 1 one can find that the spectral lags 
of the long GRBs are significantly larger than that of the 
short GRBs. The most prominent difference is observed in 

3 For a negative lag we take its minimum value as T max 
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Figure 1. Comparisons of the spectral lags between the short GRBs (solid lines) and the long GRBs: (a)the spectral lag derived from 
the whole time domain of a light curve; (b)the spectral lag of the first pulse; (c) the maximum spectral lag among pulses in a light curves. 


the t max distributions for the two kinds of GRBs. The mode 
of the T max distribution of the short GRBs is 0.015 seconds, 
while it is 0.057 seconds for the long GRBs. Soft lag is dom¬ 
inated in long GRBs. About one-third of the short bursts, 
however, have a negative spectral lag, i.e., the hard photons 
are lag behind the soft ones (hard lag). Further discussion 
on these short GRBs is referred to §6. 


5 ON THE RELATION BETWEEN THE 

SPECTRAL LAG AND THE PULSE 

DURATION 

As we shown above the spectral lags of the short GRBs 
are shorter than that of the long GRBs. One could argue 
that as short GRBs are short, all timescales will scale with 
the burst duration and the spectral lags will therefore be 
shorter. Liang et al. (2002a) showed that the pulse duration 
of the short GRBs are indeed significantly shorter than the 
long GRBs. We show the duration distributions of the pulses 
with r max for the short and long GRBs in our sample in Fig¬ 
ure 2. Their means, medians, and modes are also tabulated 
in Table 1. It is found that the pulse durations of the long 
GRBs are significantly longer and have a larger dispersion 
than that of the short GRBs. One may suspect that if there 
is a correlation between the spectral lag and the pulse du¬ 
ration. Figure 3 shows the r max as a function of the pulse 
duration for the short and long GRBs. It is found that the 
two quantities are not correlated for the short GRBs, but 
they are for the long GRBs. We measure this correlation by 
the Spearman correlation analysis method. The correlation 
coefficient is 0.76 with a chance probability p < 10~ 4 for the 
long GRBs. 

The ratio of the spectral lag to the pulse duration may 
reflect some intrinsic properties and present a clue to the 
mechanism dominated the spectral lag. We calculate this ra¬ 
tio by R = Tmax /ST. Figure 4 shows the comparison of the 
R distributions between the long and the short GRBs. It is 
found that the R distributions of the two types of GRBs 
are narrowly clustered. We examine the normality of the 



Figure 2. Distributions of the the pulse durations for the short 
(Solid line) and long (dashed line) GRBs in our sample. 

two distributions by the Shapiro-Wilk normality test. The 
normality test shows that the two distributions are normal 
in a confidence level of 3<r. We fit them by a normal dis- 

a ( fl - 0 . 023) 2 

tribution function and derive dN/dR oc e 0.076 2 and 

a ( B - 0 . 046) 2 

dN/dR <x e 0 . 076 2 for the short and long GRBs, re¬ 
spectively. The shapes of the two distributions are almost 
the same but centering at 0.023 and 0.046 for the short and 
long GRBs, respectively. 


6 THE SHORT GRBS WITH HARD LAG 

As we showed in §4, about one-third of the short GRBs 
have a hard lag. It is possible that some of them are caused 
by statistical fluctuation when the soft lags of these GRBs 
are very small. We thus consider only those short GRBs 
with T m ax < —0.01 seconds. Figure 5 shows the distri¬ 
bution of these lags. We find that the statistical fluctu¬ 
ation could not explain these hard lags. The distribution 
show a peak at ~ —0.04 seconds. It is an incomplete nor- 
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Table 1. Comparisons of the spectral lag and pulse duration distributions for the short and long GRBs in our sample 




T 


T lst 


Tin ax 

ST* 



short 

long 

short 

long 

short 

long 

short 

long 

Mean 

0.027 

0.197 

0.028 

0.209 

0.028 

0.248 

1.06 

4.01 

Median 

0.021 

0.050 

0.021 

0.056 

0.022 

0.094 

0.98 

2.00 

Mode 
Hard lag 
Pks 

0.015 

30% 

2.8 x 10“ 14 

0.023 

17% 

0.015 

31% 

2.7 x IQ" 17 

0.024 

22% 

0.015 

30% 

1.4 x 10“ 43 

0.057 

11% 

0.85 

1.25 


* The ST 

is corresponding to the 

pulse with t max in a 

burst. 






Figure 3. Correlation between the spectral lag and pulse dura¬ 
tion: (a)the short GRBs, (b) the long GRBs. The two quantities 
are not correlated for the short GRBs but they are for the long 
GRBs. The solid line in the right panel is the best fit to the two 
quantities for the long GRBs. 



Figure 5. Distribution of the spectral lags for the short GRBs 
with r m ax < —0.01 seconds. The dashed line is the best fit by a 
Gaussian model. 



R 

Figure 4. Distributions of the ratio of r max /5T for the short 
GRBs (solid step-line) and the long GRBs (dashed step-line). 
The solid and dashed lines are the fitting results by a Gaussian 
model for the short and long GRBs, respectively. 

mal distribution with a cutoff at - 0 . 01 , which can fitted by 

2 (Tma x +0.CI34) 2 

diV/(ir max <x e o.o48 2 . This likely implies that there 

are two components in the short GRB sample. The propor¬ 
tion of this component calculated by this model is ~ 17%. 
The fraction of the short GRBs with hard lags caused by 
the statistical fluctuation is thus ~ 13%, roughly consistent 


with the fraction of the GRBs with a hard lag in the long 
GRB sample. Most recently, Tanvir et al. (2005) find a po¬ 
sition correlation between the short GRBs and the nearby 
galaxies and propose a proportion of 5% ~ 25% of short- 
hard GRBs originated at low redshifts. This proportion is 
roughly consistent with that of the short GRBs with a hard 
lag. However, it is unclear if the short GRBs with a hard 
lag are indeed associated with the sources originated at low 
redshift. This is very interesting as it is probably challeng¬ 
ing for many models. We thus present a further analysis on 
these bursts. 

BATSE GRB trigger # 2365 has the maximum hard 
lag among the short GRBs. The T 90 of this burst is 1.536 
seconds. Figure 6 show its light curves in 25-55 keV and 
110 — 300 keV bands. A prominent hard photon lag is ob¬ 
served. The lag is 0.151 seconds measured by our method. 
The ratio of the rise time to the decay time of this pulse in 
25-55 keV band is 0.77 measured at the half of the amplitude 
of the light curve. This is consistent with the statistical re¬ 
sult from a sample of 92 short GRBs by Liang et al. (2002a). 
We further compare the pulse duration distribution of the 
short GRBs having a soft lag with that having a hard lag 
in Figure 7. We also could not find any peculiar features to 
distinguish both of them. 
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Figure 6. Light curves of BATSE GRB # 2365 in 25-55 keV 
(thick line) and 110-300 keV (thin line) bands. The two vertical 
lines mark the pulse with the most significant hard lag in the 
short GRB sample. 



Figure 7. Comparison of the distributions of the pulse durations 
for the short GRBs with a soft lag (thin line) and with a hard lag 
(thick line) 

7 DISCUSSION 

The precision of the spectral lag measurement is determined 
by the time resolution and the level of signal-to-noise of the 
light curves. In our analysis we derive the spectral lag by 
fitting the CCF with a Gaussian model. This method re¬ 
duces the statistical fluctuation and the time bin selection 
effect. It gives a more precise measurement of the spectral 
lag than that directly read off from the CCF, especially for 
light curves in a short time period with a lag being compa¬ 
rable or smaller than the time resolution. We use the time- 
tagged event (TTE) data to examine the reliability of our 
spectral lag measurement. The TTE data are available only 
for the short GRBs. The TTE data record the arrival time 
and energy band for each photon event. We construct the 
light curves of our short GRB sample in 32 ms time bin and 
calculate their spectral lags. It is found that both the spec¬ 
tral lags in 32 ms and 64 ms are comparable, suggesting our 
spectral lag measurement is reliable. 

The curvature effect has been extensively studies (e.g., 
Ryde & Petrosian 2002; Qin 2002; Qin et al. 2004; Ko- 


cevski et al. 2003; Dernier 2004; Zhang et al. 2006; Dyks, 
Zhang, & Fan 2005), and it is suggested to responsible for 
the spectral lag (e.g., Salmonson 2000; Ioka & Nakamura 
2001; Shen, Song, & Li 2005; Lu & Qin 2005). Ryde (2005) 
found that the spectral lag depends mainly on the pulse de¬ 
cay timescale. This is consistent with the expectation of the 
curvature effect. In this scenario, the light curve in a softer 
band tends to have a longer tail, and has a larger spectral 
lag relative to the light curve in a harder band. We thus 
intuitively expect the ratio between the spectral lag and the 
pulse duration is universal for different GRBs, if the decay 
phase of a pulse is dominated by the curvature effect. As 
we shown in Figure 3, the distributions of the ratios for the 
short and long GRBs are normal with the same width, sug¬ 
gesting the ratios are universal within the long GRBs and 
short GRBs, respectively. However, they peak at different 
values. The narrowly clustering of the distributions for the 
two types of GRBs hints that the curvature effect may be 
partially responsible for the observed spectral lag. The shift 
of the peaks of the two distributions, however, implies that 
this effect alone could not account for the difference of the 
spectral lags between the short and long GRBs. Other effects 
should play an important role on the pulse evolution, making 
the difference of the spectral lags in the two types of GRBs. 
Several scenarios are also involved to present an explanation 
for the spectral lag, such as an intrinsic cooling effect of the 
radiating electrons (e.g. Zhang et al. 2002) and a Compton 
reflection of a medium (e.g., Ryde & Petrosian 2002). An¬ 
other possible scenario might be related to the activity of 
the central engine and hydrodynamic timescale of the in¬ 
ternal shocks (Daigne & Mochkovitch 1998; 2003). If the 
shell thickness is significantly smaller than the initial shell 
separation, the hydrodynamic timescale may be relatively 
short enough and the curvature effect governs the pulse evo¬ 
lution. If the shell thickness is comparable to the initial shell 
separation or the outflows are continuous, the hydrodynam- 
ical timescale may dominate the pulse evolution. Daigne & 
Mochkovitch (1998; 2003) developed a model in the frame¬ 
work of internal shock model (Rees & Meszaros 1994) and 
found the observed relation between luminosity and spectral 
lag (Norris et al. 2000) indeed can be reproduced when the 
hydrodynamical timescale is taken into consideration. 


The short GRBs with a hard lag is very interesting. This 
phenomenon is a challenge to many models. The proportion 
of the short GRBs with negative lag is ~ 17%, excluding the 
statistical fluctuation effect. We notice a new finding by Tan- 
vir et al. (2005), who found a position correlation between 
the short GRBs and the nearby galaxies and proposed a pro¬ 
portion of 5 ~ 25%. This proportion is roughly consistent 
with that of the short GRBs with a negative lag. However, it 
is unclear if the short GRBs with negative lag are indeed as¬ 
sociated with the sources originated at low redshift. Current 
GRB models are hard to explain this phenomenon. We have 
tried to identify some distinguished features of these bursts 
from those with a soft lag. We do not find such signatures 
from their duration and the ratio of the rise time to the 
decay time. We guess that a more physical mechanism dom¬ 
inated the hard lag may be hid behind the morphological 
features of the light curves. 
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8 CONCLUSIONS 

We have present a detail analysis on the spectral lags of 
the short GRBs and compare them with that of the long 
GRBs by using the CGRO/BATSE GRB Catalog. Our sam¬ 
ple includes 308 short GRBs and 1008 long GRBs. The light 
curves of these GRBs are in 64 ms time bin and they have at 
least one long and intense pulse, which satisfies 5T > 0.512 
seconds at c = la and c max > 6cr. We calculate the CCFs 
for the light curves in 25-55 keV and 110-300 keV bands and 
derive the spectral lags by fitting the CCF with a Gaussian 
model. We summary our results as follows. 

(1) The spectral lag distribution of the short GRBs is 
significantly different from that of the long GRBs. The null 
hypothesis that the spectral lag distributions of the two 
kinds of GRBs are from the same parent is rejected, which 
is tested by the K-S test. 

(2) Excluding the statistical fluctuation effect, a propor¬ 
tion of ~ 17% of the short GRBs have a negative spectral 
lag. We do not find any peculiar features from their light 
curves to distinguish these bursts from those with a posi¬ 
tive spectral lag. We argue that a more physical mechanism 
dominated the hard lag may be hid behind the morpho¬ 
logical features of the light curves. This should be a great 
challenge to the current GRB models. We notice that this 
proportion is consistent with the proportion of short GRBs 
correlated with nearby galaxies newly discovered by Tanvir 
et al., although it is unclear if these short GRBs are indeed 
associated with the sources originated at low redshift. 

(3) While the spectral lags of the long GRBs are strongly 
correlated with the pulse duration, the spectral lags of the 
short GRBs are not. However, the distributions of the ratio 
Tmax/ST for two kinds of GRBs are normal in a confidence 
level of 3a, and they have the same shapes but centering at 
0.023 and 0.046 for the short and long GRBs, respectively. 
This result indicates that the curvature effect alone could 
not explain the difference of the spectral lags between two 
types of GRBs. The radiative processes at work in GRBs and 
the hydrodynamic timescales of the outflows may also play 
an important role as suggested by Daigne and Mochkovitch 
(1998, 2003). 
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tion of China ( No.10463001 and No.10273019) and the Spe¬ 
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